function diff = model_SA( param, aux)
    
    c = fzero(@ (x) model_SA_inner( x, param, aux), 1, optimset('display','off'));
	param.c = abs(c);
    diff = model_2( param, aux);
    
end